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J&1AGNETIC RESONANCE Jjg?!HQD_WI.TH^ MON- LINEAR MAGNETIC FIELD GRADIENTS 



The invention relates to a magnetic resonance (MR) method for the imaging of 
an object arranged in a steady magnetic field, whereas the following steps being executed 
according to said method: 

- applying of a stationary magnetic field and temporary magnetic position-dependent field 
5 patterns, 

- excitation of spins in a part of the object, 

- acquiring magnetic resonance signals by one or more receiver antennae, 

- and MR signals are acquired during application of the position-dependent field patterns. 

The invention also relates to an MR device for carrying out such a method. 
10 Current MR systems are usually very limited in size of the useable imaging 

region which is dependent on the homogeneity volume of the main stationary magnetic field. 

In this volume the main magnetic field has to obey stringent aims on constancy and the three 

temporary magnetic gradient fields have to obey stringent aims with respect to linearity. 

Increasing the volume of the main magnetic field, where the field is substantially constant 
15 and the normal gradients are substantially linear, makes the system extremely expensive. 

It is thus an object of the present invention to obtain an MR imaging system 

which is less dependent on the homogeneity of the main magnetic field and of the linearity of 

the gradient fields. 

This object of the invention are achieved by a method as defined in Claim 1. 
20 The invention is further related to an apparatus as defined in Claim 6 and to a computer 
program product as defined in Claim 9. 

The main advantage of the present invention is that the method allows the use 
of any position-dependent field pattern, which is substantially non-linear, so that subsampiing 
method like SENSE can also be applied for a MR systems with imperfect coils. 
25 These and other advantages of the invention are disclosed in the dependent 

claims and in the following description in which an exemplified embodiment of the invention 
is described with respect to the accompanying drawings. Therein shows: 
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Fig. 1 diagrammatically a magnetic resonance imaging system in which the 
present invention is used, 

Fig. 2 a one-dimensional object acquired with two field-patterns, 
Fig. 3 the x-component of the field caused by the two field patterns, 
5 Fig. 4 the locus of the object-space in G-space, 

Fig. 5 folding of the object-locus in G-space due to undersampling, 
Fig. 6 the position dependent field patterns of four conductor systems, 
Fig. 7 the locus of the object in G-space, 

Fig. 8 the same locus as in Fig. 7 with a single measurement point, and 
10 Fig. 9 the same locus as in Fig. 7 with several measurement points. 



Fig. 1 shows diagrammatically a magnetic resonance imaging system in which 
the invention is used. The magnetic resonance imaging system includes a set of main coils 10 
1 5 whereby the steady, uniform magnetic field B 0 is generated. The main coils are constructed, 
for example in such a maimer that they enclose a tunnel-shaped examination space. The 
patient to be examined is slid into this tunnel-shaped examination space. The magnetic 
resonance imaging system also includes a number of field pattern coils 11,12 whereby 
magnetic fields exhibiting spatial variations, notably in the form of temporary position- 
20 dependent field patterns in individual directions, are generated so as to be superimposed on 
the uniform magnetic field. The field pattern coils 1 1, 12 are connected to a controllable 
power supply 21. The field pattern coils 11, 12 are energized by application of an electric 
current by means of the power supply unit 21. The strength, direction and duration of the 
field patterns are controlled by control of the power supply unit. The magnetic resonance 
25 imaging system also includes transmission and receiving coils 13, 15 for generating the RF 
excitation pulses and for picking up the magnetic resonance signals, respectively. The 
transmission coil 13 is preferably constructed as a body coil whereby (a part of) the object to 
be examined can be enclosed. The body coil is usually arranged in the magnetic resonance 
imaging system in such a manner that the patient 30 to be examined, being arranged in the 
30 magnetic resonance imaging system, is enclosed by the body coil 13. The body coil 13 acts as 
a transmission aerial for the transmission of RF excitation pulses and RF refocusing pulses. 
Preferably, the body coil 13 involves a spatially uniform intensity distribution of the 
transmitted RF pulses. The receiving coils 15 are preferably surface coils 15 which are 
arranged on or near the body of the patient 30 to be examined. Such surface coils 15 have a 
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high sensitivity for the reception of magnetic resonance signals which is also spatially 
inhomogeneous. This means that individual surface coils 15 are mainly sensitive for 
magnetic resonance signals originating from separate directions, i.e. from separate parts in 
space of the body of the patient to be examined. The coil sensitivity profile represents the 
spatial sensitivity of the set of surface coils. The receiving coils, notably surface coils, are 
connected to a demodulator 24 and the received magnetic resonance signals (MS) are 
demodulated by means of the demodulator 24. The demodulated magnetic resonance signals 
(DMS) are applied to a reconstruction unit 25. The reconstruction unit reconstructs the 
magnetic resonance image from the demodulated magnetic resonance signals (DMS), on the 
basis of the coil sensitivity profile and the knowledge of the fields of the temporary field 
patterns of the set of surface coils. The coil sensitivity profile is stored, for example 
electronically, in a memory unit which is included in the reconstruction unit. The 
reconstruction unit derives one or more image signals from the demodulated magnetic 
resonance signals (DMS), which image signals represent one or more, possibly successive 
magnetic resonance images. The reconstruction unit 25 in practice is preferably constructed 
as a digital image processing unit 25 which is programmed so as to reconstruct the magnetic 
resonance image from the demodulated magnetic resonance signals and on the basis of the 
coil sensitivity profiles and the knowledge of the fields of the temporary field patterns. The 
image signal from the reconstruction unit is applied to a monitor 26 so that the monitor can 
display the image information of the magnetic resonance image (images). It is also possible 
to store the image signal in a buffer unit 27 while awaiting further processing, for example 
printing in the form of a hard copy. 

In order to form a magnetic resonance image or a series of successive 
magnetic resonance images of the patient to be examined, the body of the patient is exposed 
to the magnetic field prevailing in the examination space. The steady, uniform magnetic field, 
i.e. the main field, orients a small excess number of the spins in the body of the patient to be 
examined in the direction of the main field. This generates a (small) net macroscopic 
magnetization in the body. These spins are, for example nuclear spins such as of the 
hydrogen nuclei (protons), but electron spins may also be concerned. The magnetization is 
locally influenced by application of the position-dependent field patterns. Subsequently, the 
transmission coils apply the RF excitation pulse to the examination space in which the part to 
be imaged of the patient to be examined is situated. The RF excitation pulse excites the spins 
in the region of interest, i.e. the net magnetization then performs a precessional motion about 
the direction of the main field. During this operation those spins are excited which have a 
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Larmor frequency within the frequency band of the RF excitation pulse in the main field. 
However, it is also very well possible to excite the spins in a part of the body which part is 
much larger than a thin slice; for example, the spins can be excited in a three-dimensional 
part which extends substantially in three directions in the body. After the RF excitation, the 
spins slowly return to their initial state and the macroscopic magnetization returns to its 
(thermal) state of equilibrium. The relaxing spins then emit magnetic resonance signals. 

In contradistinction to a conventional MR imaging system as described in 
principle above the MR imaging system according to the present invention contains N 
subsystems that can induce purposely and substantially position-dependent field patterns. The 
expression "field pattern" is meant to be in more general terms for what is conventionally 
mentioned "gradient" or "main field" which are normally linear and homogeneous. The field 
patterns used in this invention are substantially non-linear. A further important feature of the 
method is that the number N of field patterns is always larger than 3. In addition, at least N - 
1 field patterns are independently controllable in their field strength, i.e. the strength of the 
created pattern is controlled by the current applied to the particular subsystem, which is 
called here "switchable". 

In the sequel, these field patterns are called G 0 (x, y, z), Gi(x, y, z), G 2 (x, y, z), 
Gn_i(x, y, z), or for simplicity G 0 , Gi, G 2 , Gn-i- All of them are switchable, except the 
main field pattern Go which is comparable with the main field Bo in a conventional MR 
system. A conventional MR system may be regarded as one where Gi(x, y, z) = G x (x, y, z) = 
x, ... , G 3 (x, y, z) = G 2 (x, y, z) = z and there is no other purposely position-dependent field. 

The acquisition of data on such a system is performed by using any known 
MR method as described above in connection with Fig. 1, whereas Gi, G 2 , G3, G4, .... are 
switched as in "normal" MR acquisition systems for G*, G y and G z . As a main difference, this 
is an experiment of higher dimensionality, e.g. using four encoding dimensions to acquire 
information from a three-dimensional human body. In order to compensate for this gross 
increase in scanning time as more data must be sampled, a further important characteristic of 
the present method is that the acquired data is grossly undersampled in analogy to the well- 
known SENSE technique (see e.g. K. Pruessmann et. al. in Proc. ISMRM, 1998, abstracts pp. 
579, 799, 803 and 2087). 

The reconstruction is formed by following elements: 
1. The knowledge of the mapping of co-ordinates (x, y, z) in real space to G-space. For 
better understanding of this notion suppose that each of the subsystems is driven by a unit 
current (e.g. 1A), and, if applicable the steady-state subsystem G 0 is driven by its standard 
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current. Then every point in real space will experience a given field component caused by 
the system Go, some other field component caused by the system Gi, etc. Thus any 
specific point (x, y, z) will be mapped to a point (go, gi, g2, gN-i) in G-space, as is 
depicted in Fig. 4. The mapping is, in principle, known from the system design. 

2. The acquired raw MR-data can be seen as an N-dimensional set of data in the (ko, ki, k2, 

k N .0-space in analogy to to (k x , k y , k z ). That data can be transformed by Fourier 
transformation to G-space (cf. "image space"). 

3. Some pairs of (x, y, z) locations wiU map to identical (go, gi, g2, .» gN-i) - co-ordinates 
due to folding. These ambiguities can be resolved by knowledge of the coil-sensitivity 
profiles which technique is known from SENSE, or by continuity constraints or by known 
sparsity of the object, e.g. when vessels are imaged. 

For better understanding of the present method consider now the simplified 
case of an MR imaging system designed to image an object in only one dimension. The 
imaging system is depicted in Fig. 2 wherein Go is a circular coil system and Gi is a straight 
lead that is close to the one-dimensional object 31, i.e. (part of) the patient. It is further 
assumed that a strong constant field is present, aligned with the main direction of the object 
31. In Fig. 3 the field strength B of both field patterns G 0 and Gi are depicted with respect to 
the x-axis. Given these field patterns, one can draw the locus 32 of the one-dimensional 
object in the two-dimensional G-space as shown in Fig. 4. 

In principle on this system a normal 2-dimensional MR experiment can be 
applied, whereas the current in Go will be driven as a readout field pattern and the current in 
Gi as a phase-encoding field pattern. After Fourier transformation a 2-dimensional image 
will be obtained, of which only a small part would be nonzero, i.e. the locus of the object 31. 
However, such a simple approach would have the significant drawback of a gross increase in 
scanning time. In above example 256 2 measurements points are needed to resolve 
approximately 256 points of the patient. In order to overcome this disadvantage, a gross 
subsampling of G-space can be applied. The result can nevertheless be reconstructed, as large 
areas of G-space are factually empty. In Fig. 5 the folding of the locus 32a, 32b, 32f of 
object 31 in G-space due to undersampling is shown. Large parts of the object 31 can be 
uniquely extracted from G-space, but on some points there may be an overlap. This can be 
intrinsic as shown at circle 33 or caused by undersampling as shown at full bullet 34. The 
intrinsic ambiguity may be caused by system design (e.g. due to current lead constraints, two 
points of the object have exactly the same magnetic field under all circumstances etc.). 
However, it is more likely that the ambiguity is caused by undersampling. Both ambiguities 
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can be resolved in several manners. The first one is to apply the knowledge on coil sensitivity 
patterns, if the object has been acquired using an array of coils - like in the SENSE method 
mentioned above. A second solutions is the use of continuity constraints e.g. by assuming 
that is unlikely that an object will show a distinct local detail at exactly a particular cross- 
over position. Especially if such a system is used for angiographic purposes, the knowledge 
can be used that large parts of the interesting object will contain only very little signal at all. 

This more theoretical approach can be understood better by way of the 
following model as depicted in Fig. 6. If a more or less linear field pattern system in direction 
"feet-head" (z-direction) would be used like e.g. the open magnet system as shown and 
described in Fig. 4 of WO-A-00/33100, the field pattern system must store a huge amount of 
energy which only can be arranged by very expensive electric amplifiers, which imply very 
strong field variations at about 1 meter off the centre of the MR system. Therefore, nervous 
stimulation will occur with the patient and the operator. It should be clear that such large 
fields with high energy storage will make the practical use of linear field patterns in z- 
direction almost impossible. An alternative to a linear field pattern in z-direction can be 
obtained by a system of parallel conductor systems, e.g. the four systems A, B, C and D, each 
of which has a design similar as the main magnet field of the system as shown in Fig. 2 of 
WO-A-00/33 100, whereas the "switched" magnetic fields are substantially shorter in z- 
direction. The four conductor systems A, B, C and D are identical, but provided equidistantly 
at different positions along the z-axis. Each of the conductor systems A, B, C and D can be 
controlled separately. The field pattern of each of the conductor systems is shown in Fig. 6, 
whereas Bgs is the extra field "switched" by the conductor system at a unit of current. Since 
the field pattern of the conductor systems is a Gaussian function and thus not linear in z- 
direction, it is definitely not a "gradient" as known from the usual MR systems. However, the 
conductor systems can be used similar as "gradients" in the classical manner since they have 
a pattern which varies continuously in z-direction and can be switched on and off. It is now 
possible to perform two different MR sequences in z-direction, whereas in x- and y-direction 
normal gradients are used. The first sequence the combination Gi = (A + B - C - D) is 
applied and the second sequence the combination G2 = (- A + B + C — D) is applied. For each 
of the sequence the data received by the receiver antennae will be reconstructed in a normal 
maimer, e.g. according to a 3-dimensional Fast Fourier Transformation (FFT). For each of 
the (x, y) coordinate of the patient a graph in Gi - G 2 space according to Fig. 7 will be 
obtained which is commensurate to a straight line parallel to the z-axis. Concerning the axis 
Gi and G 2 one should realise that it represents the magnetic field of the combination of coils 
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or conductor systems A, B, C and D excited by a unit of electric current. So, the first MR 
sequence is measuring the equidistant points along this axis (for all x and y coordinates of 
that point along the z-axis). The curve in Fig. 7 is the locus of all combinations which are 
possible within the patient's body. For instance, the point at the far left of the patient (circle 

5 35 in Fig. 6) which is almost only excited by conductor system A, experiences a negative 
magnetic field for G 2 and a positive magnetic field for Gi. The results of the mentioned 
measuring sequences represent in fact a projection of the patient curve on the axis. Thus, the 
measurement of point 36 ("X") at the Gi axis is a sum of images within the patient at two 
different positiions 37 ("full bullets"), as can be seen in Fig. 8. The two points can be derived 

1 0 from the measurement points 36 by knowledge of the sensitivity of the RF coils in analogy 
to the SENSE method. Since the number of measured points on the curve is essentially 
identical to the number of unknown variables, the calculation should be performed with all 
other measured variables. In the example showed in Fig. 7 the number of measurement points 
37 ("full bullets") is seven and also the number of unknown variables 36 ("X") is seven - see 

15 Fig. 9. 

With the system described above one can design an MR system with imperfect 
coils having highly non-linear field patterns - comparable with "gradients" in classical terms 
- in order to obtain a relatively inexpensive system of large field-of-view. Important 
prerequisites for such a system are that the position-dependent field patterns and the coil 
20 sensitivity patterns, i.e. the proper features of the RF coil system, should be known. 



